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, ^ ' Abstract 

I Using the covariant spectator theory and the transversity formahsm, the unpolarized, coincidence 
cross section for deuteron electrodisintegration, d{e, e'p)n, is studied. The relativistic kinematics 

^ . are reviewed, and simple theoretical formulae for the relativistic impulse approximation (RIA) are 

OO , derived and discussed. Numerical predictions for the scattering in the high region obtained 

I from the RIA and five other approximations are presented and compared. We conclude that 

Tij" ' measurements of the unpolarized coincidence cross section and the asymmetry A^, to an accuracy 

■ that will distinguish between different theoretical models, is feasible over most of the wide kinematic 

O ■ range accessible at Jefferson Lab. 

PACS numbers: 21.45.+v,25.10.+s,25.30.Fj 

> 
X 

C3 ; I. INTRODUCTION 

Inelastic scattering of electrons from the deuteron is an important source of information 
about the nuclear current, deuteron structure, and the A^A^ force. The exclusive scattering 
cross section, d{e,e'p)n, was first measured almost forty years ago and since then it 
has been measured under a wide variety of kinematic conditions p|. There is a substantial 
body of data for this reaction, including cross section measurements [0, |, |^, || as well as 
separations of various response functions 0, ||, ^ |10|, [11], |12], |13|, |14|, |15|, |16| which differentiate 
between absorption of longitudinal and transverse photons. 

In this paper we survey results that might be expected from a new generation of d{e, e'p)n 
coincidence measurements proposed for Jefferson Laboratory (JLab). At JLab it is possible 
to carry out a comprehensive program of measurements at both high and large W (where 
W is the invariant mass of the final np state). A broad program of such measurements offers 
the best hope of independently determining effects of final state interactions and the nuclear 
current, permitting the extraction of important new information about the short range A^A^ 
interaction. 
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Electro disintegration of the deuteron has been studied theoretically by many groups. 
Recently, Arenhovel, Beck, and Wilbois [l7| have emphasized that the relativistic effects in 
inelastic scattering can be very large, even at modest momenta, and it is therefore partic- 
ularly important to have a fully relativistic theory available for the analysis of the higher 
momenta data that will be measured at JLab. Relativistic calculations of this reaction date 
back to the early work of Durand JTB| and McGee and lead up to more recent work by 



Tjon PU|. One of the goals of this paper is to present a fully modern, covariant treatment 
of this process suitable for the analysis of JLab data. 

This paper imbeds the dynamical calculation in the general formalism developed in 



Ref. [^, where a covariant, systematic treatment of most of the polarization observables that 
can be measured in the d{e, e'p)n reaction were classified and defined. There it was found 
that the use of transversity amplitudes (closely related to helicity amplitudes) gave a very 
efficient description of the reaction. Transversity amplitudes were discussed by Moravcsik 



22 1, who found that they build in the constraints imposed by parity and rotational invari- 



ance in the most efficient way. In a transversity basis, the constraints imposed by these 
symmetries insure that half of the possible amplitudes vanish identically, so that the cum- 
bersome linear relations needed in other formalisms ||23| are unnecessary. This economy will 
be essential some time in the future when large data sets exist, and it may be important to 
know whether or not a proposed new measurement will really be independent of amplitudes 
already measured. 

The details of the calculation are carried out using the covariant spectator theory, which 
has been successfully applied to the description of NN scattering and the electromag- 
netic form factors of the deuteron |^ . One feature of this theory is that the deuteron bound 
state is described by the covariant dnp vertex with one nucleon on mass-shell, and this is 
precisely the amplitude that is needed for the relativistic impulse approximation (RIA), 
making the theory well suited to the analysis of the d{e,e'p)n coincidence reaction [Q. In 
this first application of the covariant spectator theory using the transversity basis, we present 
the RIA calculation only. This provides the opportunity to work out several new technical 
details for the simplest case, and to compare to other approximations. The inclusion of final 
state interactions and interaction currents will be the subject of future work. 

A second purpose of this paper is to estimate the size of the unpolarized d{e, e'p)n cross 
sections expected over the broad range of and W accessible to JLab. In preparing 
this survey we found that relativistic and nonrelativistic predictions for d{e, e'p)n at high 
Q^, where the cross section is most sensitive to the theory, often differ by as much as an 



order of magnitude, extending the observations of Arenhovel, Beck, and Wilbois |T^. Since 



nonrelativistic calculations cannot be taken seriously at such high energies, we report our 
results for a variety of relativistic or semi-relativistic models only. In this first exploratory 
study the goal is to provide only a rough survey of the landscape. The simplicity of the 
RIA allows a uniform treatment over the entire kinematic range, but is, of course, very 
incomplete. It is our intention to follow up this study with complete calculations for cases 
where the theoretical effects look especially interesting. 

Our notation for the cross section and the spectator theory for the RIA are reviewed 
briefly in Sec. II, numerical results presented in Sec. Ill, and conclusions given in Sec. IV. 
Many theoretical details are given in the several Appendices, which are an important part 
of this paper. 
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FIG. 1: The kinematics of electron scattering when the final hadronic state is broken into two fragments 
with momenta pi and p2- 

II. THEORY 

In this section we define the coincidence cross section and the RIA matrix element. All 
other theoretical details can be found in the Appendices. 

A. The Cross Section 

Figure |1] shows the kinematics for the process e + d^e'+p + n (using the notation of 
Ref. |2^). The incident and scattered electron momenta form a plane called the "scattering 
plane" while the momenta of the proton and neutron in the final state form a second plane 
called the "ejectile plane." The virtual photon momentum is common to the two planes 
and is chosen as the direction of the 2;-axis. The two planes, which are represented by the 
{x, y, z) and (x', y', z') coordinate systems, are oriented at a relative azimuthal angle of (j). A 
rotation of the response tensor (defined below) from the unprimed to the primed frame can 
be used to extract all of the dependence from the tensor. Using this along with the explicit 
form of the electron tensor, the cross section can be shown to be of the form (c.f. Eq. (95) 



of Ref. H) 




r'-P + stB!P - - [cos 20^?}. + sin 20^^ 



dQldE'dY. 



+slt cos (puij- + sin (pRiji + 2h st'Rt' 




(2.1) 



where 



a cos \6 



(2.2) 



2E single 
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is the Mott cross section and dH is defined below. The quantities 6, E, E', and fl' are the 
electron scattering angle, the energies of the initial and final electron, and the solid angle 
of the final electron, all in the lab frame. The deuteron mass is Md and h = ±1/2 is the 
helicity of the incident electron. The electron kinematical factors are 

ST ='^+e s,T= -^(i+e^' 



V2 



where 



ST' =e(l + a^ SLT'= (2.3) 

^ = %tan- (2.4) 
with 

2 ,n2 2 2 2 2 /o 

q =-Q =v - = z/q - go (2.5) 

the square of the virtual photon four-momentum, with {z^, q^} and {t'cQo} the energy 
and three-momentum of the virtual photon in the lab and cm. systems, respectively, and 
go = |qo|, etc. 

There are two inertial reference frames that are of interest in the calculation of deuteron 
electrodisintegration: the laboratory frame which coincides with the rest frame of the target 
deuteron, and the center of momentum (cm.) frame in which the total three-momentum of 
the final state proton-neutron pair (or of the initial virtual photon and the target deuteron) 
is zero. One of the virtues of Eq. ( |2.1| ) is that the response functions R are covariant, and 
hence ( ^.Ij ) can be used to describe the cross section in either the cm. of the outgoing np 
pair or the laboratory frame by the replacement of dT, by 

dJ:\cm = P dn* (2.6) 

in the cm. frame or 

dJ:\iab = Pi dQi n (2.7) 

in the laboratory frame. [Except for special notation used in Eq. ( |2.5| ), we use a roman 
character for the magnitude of a three- momentum, so that Pi=|pi|, to distinguish it from 
the corresponding four- momentum, pi] The factor 

W 1 

^ = T:r -/ F TT-r- (2-8) 

Md . , t/pi -Eigcos^i X ^ ' 

I + MaVi )l 

is the recoil factor, where W is the invariant mass of the outgoing pair and the subscript L 
means that each variable in the parentheses is to be replaced by its value in the lab frame 
(for example, g g^ and pi is a function of 6i, the angle between the outgoing proton and 
the z axis in the lab system). We will sometimes use an asterisk (*) to denote a variable 
in the cm. system. Notation for some of the most important variables is summarized in 
Table |. 

The cm. frame (referred to as the "antilab" frame in Ref. [0) is of interest for theo- 
retical reasons because integrating over the final state kinematical variables is particularly 
convenient in this frame and the partial wave expansion of the final state is normally carried 
out in this frame. While this partial wave expansion is particularly convenient at low and 
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TABLE I: Notation for frequently used variables. 



variable 


Lab 


cm. 


photon energy 


V 




magnitude of photon 3-momentum 


Ql 


Qo 


deuteron 4-momentum 


P 


p* 


deuteron energy 


Md 


Do 


proton 4-momentum 


Pi 


Pi 


neutron 4-momentum 


P2 


P*2 


proton angle 


Ox 


9* 


neutron angle 


02 


e* + TT 


magnitude of proton 3-momentum 


Pi 


p 


magnitude of neutron 3-momentum 


P2 


p 



medium energies of a few hundred MeV, we would like to point out that the partial wave 
approach becomes extremely tedious and/or impractical at GeV energies. At such energies 
Glauber theory |^^, or the new so-called "three-dimensional" methods of calculating the 
NN amplitude directly without partial wave expansions |2^, are better. In any case, since 
the final scattering state is by far the most complicated ingredient in the calculation of the 
transition matrix elements, it is important to be able to carry out calculations in this frame 
and to translate them to the lab frame. The necessity of boosting the calculation from the 
cm. frame to the lab frame requires that the Lorentz properties of the matrix elements be 
understood. This goal is conveniently accomplished by using the Jacob and Wick helicity 
formalism provided that it can be shown that the various ingredients in the calculation 
of the matrix elements, such as the wave functions, are covariant. We will assume for the 
moment that this is the case and will demonstrate later that it is true for the particular 
calculations which are described in this paper. 

The nine response functions of (|2.1| ) are related to sums over the squares of matrix 
elements of the deuteron current. In the helicity basis, with the helicity of the virtual 
photon, Ai and A2 the helicities of particles 1 and 2 in the final state, and A^ the helicity 



-^A^(^) ^d)- Following the 



of the initial deuteron, the current operator is written <^AiA2 
conventions of Jacob and Wick we choose particle 1 in the 
and particle 2 to be the neutron. The current operator conserves parity, which means that 
the matrix elements satisfy the condition 



:inal state to be the proton 



(A1A2 \JxM\ ^d) = ± (-Ai, -A2 \J-xM\ -^d) 



(2.9) 



where the phase depends on the helicities (see Ref |2T]]). For this reason it is convenient to 
introduce symmetric and antisymmetric combinations of the |A^| = 1 amplitudes 

"/'a, A,(Pi,P2, q) = (A1A2 \Js{q)\ \d) = \{{\i\2 \Mq)\ \d) - (A1A2 | J-i(g)| A,) } 

^a\A,(Pi,P2,g) = (AiA2|Ja(g)|A,) = ^{(AiA2|Ji(g)|A,) + (AiA2|J-i(g)|A,)}, (2.10) 

where, because of the phases, Jg is symmetric under the Y parity transformation (parity 
followed by rotation by n about the y axis) and Ja is antisymmetric (and we note for future 
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TABLE II: Response functions. 



R 



00 



(I) 

L 

R^p = 2{Raa + Rss) 
Rj'j' = '^{Raa Rss) 



(I) 
LT 



R 

o(I) 



4Re Ros 
-- Aim Ros 



Rm ^ — 4Re Rs 



(11) 

TT 
(11) 
LT 

5(11) 
LT' 



R. 
R 
R 



-4Im Rsa 



4Im R 



Oa 



4Re Roa 



reference that Jqx^x,^ = (A1A2 |^o(q')| ^d) is also symmetric). We then define the deuteron 
response tensors R^gg, and i?^"? 



,2 



.2 



^S' = ^ E Ej(p'^)A;A.J,\A,(pi,P2,g)(pi,-''^)A,A:,4'4,,(pi,P2,^ (2.11) 

where 5^ and (7' = {0, s, a}, and and p^) are the spin density matrices for one nucleon in the 
final state or the deuteron target, with p"*" being the part of the density matrix symmetric 
under Y parity and p~ the part antisymmetric under Y parity. Symmetry under the Y 
parity operation then insures that those observables of type (II) must include one, and only 
one factor of the antisymmetric current Ja (further details can be found in Ref. |2l|]). The 
relation between the nine response functions that appear in Eq. (p.l|) and the tensors defined 
in (|2.11|) are given in Table |I[ The normalization of Eq. ( p.ll|) and the density matrices is 



consistent, for unpolarized reactions, to summing over final state spins and averaging over 
the initial deuteron spin. 
For unpolarized particles, 

(P7v)a;ai = ^5a;ai (p^)a;ai = 

ipt)x,K = l^y, (pB)a,a— 0, (2.12) 

so the observables of type (II) are zero. If we also limit discussion to unpolarized electrons, 
the terms proportional to the electron helicity h average to zero, and the cross section 
depends on only 4 response functions: 



1^?^ + stR^ - ^ cos 20 R^T + slt cos ^?^| (2.13) 



dn'dE'dL AnMd ql 

A second, independent combination of the same four response functions gives the asymmetry 

d^a . , (i^cr 
4 _ dfl'dE'dE ^'^ ^ ~ dQ'dE'di: = _ sltR^t 

't^ (0 = 0) + — U = vr) + StRt - 2^TT 

dn'dE'dj: ^ dn'dE'dj: ^ 
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FIG. 2: Left panel is the Q'^v plane, and the right panel the W^v plane. In each panel lines with constant 
X are shown, and the cases analyzed in the following section are shown as dots. The shaded area in each 
panel is the region where pion production is kinematically possible. Note that pions cannot be produced 
near the line a; = 2, but that inelasticity sets in even at small along the quasielastic ridge (a; = 1) and 
at smaller values of x. The dashed line in the left panel corresponds to W'^ = 9m^, just above the region 
shown in the right panel. 



where the electron kinematics is held fixed and the outgoing proton is measured forward to 
the direction of the virtual photon momentum q (at = 0) and backward (at = vr). The 
longitudinal contributions e!j} — i^Rj)j. can be separated from the transverse response Rj^ 
by measuring the cross section for the same kinematics at forward and backward electron 
scattering angles, but the transverse interference term i?^^ can be separated from e!^' only 
by an out-of plane measurement (for example, (p = tt/2). 

These four unpolarized structure functions are only a small fraction of the structure 
functions which can be measured. With polarized electrons, targets, and recoiling nucleons 
many more can be studied 1^1], but these observables tend to be very sensitive to final 
state interactions and interaction currents. In this first paper we have omitted final state 
interactions and interaction currents, and hence also limit the discussion to unpolarized 
observables. 



B. Kinematics 

The response functions R depend on three variables: Q^, and the angle 9i between pi 
and q, where pi is the three-momentum of the particle detected in coincidence with the final 
electron (assumed here to be the proton). The variables and v are fixed by the virtual 
photon, and we choose 9* (the lab value of ^i) rather than 9i because it is independent of 

and v and always varies between and vr. In place of z/, it is often convenient to use 
oi X = Q"^ /2mv, the Bjorken scaling variable. The mass of the final state, W , is related to 
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FIG. 3: The relativistic impulse approximation (RIA) to deuteron electrodisintegration in the lab frame 
and the cm. frame. The open circle denotes the deuteron wave function, the filled circle the nucleon form 
factor, and the nucleon propagating between the two is off-shell. Note that the wave functions ip always 
have one particle off-shell. 



z/ (or x) by 

= M',+2M,u[l-—J (2.15) 

The region of allowed values of and u is shown in Fig. |^. If the scattering is elastic, so 
that the deuteron remains bound after the scattering, x ~ 2 , and this defines one boundary 
of the allowed scattering region. It is sometimes assumed that pions must necessarily be 
produced as increases, but as long as x remains close to 2, the final state remains below 
the pion production threshold up to very large values of Q^, and one may try to explain the 
large behavior of these inelastic processes using a theory with no pion rescattering in 
the final state. The line x = 1 is the quasielastic peak; when is large the region between 
X = 1 and X = is the region where y (or x) scaling is observed. If x is small, pions will 
be produced more and more easily as increases (penetrating further and further into 
the shaded region in Fig. and explicit treatment of the pion degrees of freedom will be 
necessary. 

The variables Q^, v (or x or and Q* are convenient for thinking about final state 

interactions. However the RIA depends primarily on only two variables, and Pmiss, where 
Pmiss is the value of the spectator momentum in the lab system. The spectator momentum 
Pmiss may be either pi or p2 depending on which of the two nucleons was struck by the virtual 
photon (in the absence of final state interactions or interaction currents, this is all that can 
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FIG. 4: Polar plots showing the locus of the momentum vector pi in the lab system (solid lines are 
relativistic; dashed lines nonrclativistic) . The horizontal axis in each panel is z; the vertical is x. The three 
left hand panels have = 1 GeV^ and the various values of x and q shown on each panel; the right hand 
panels are for = 3 GeV^. In all panels the q vector points to the right along the z axis and sets the 
scale. The other two vectors are pi and P2 for the symmetry case discussed in the text. 
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happen), and the cross section is therefore the coherent sum of two terms. Symbohcally, the 
RIA current is 2 

jRIA(pi,P2,g) = |F^(gXp2)±F"(g2)V^(pi)| , (2.16) 

where exphcit formulae for the magnitudes of the rest frame three momenta, pi and p2, are 
given in Eq. ( |2.19| ) below. The two terms contributing to this sum are illustrated in Fig. |^. 
Since is normally a rapidly decreasing function of p, these two terms are normally 

dominated by the one with the smallest Pmiss- 

The momenta pi and p2 are most easily obtained by boosting from the cm. frame. Their 
X and z components are 

= ±p sin 9* 

P-= |±P^cosr (2.17) 

where the upper (lower) sign is for i = 1 (i = 2), p is the magnitude of the nucleon momenta 
in the cm. frame, Ew is the energy of the outgoing pair in the lab frame, with 



P 



2 



Ew= ^W^ + ql=Ma + u, (2.18) 
and the other variables were previously defined (recall Table |I|) . Hence 

P.^=(|±p|^cosr)' + p2sin2r (2.19) 

The behavior of the magnitudes of pi and p2, and the angles 61 and 62, for six choices of 
and x, can be inferred from Fig. |^. The solid lines in each panel are the locus of points 

swept out by Eq. ( [^.17D , and the dashed lines by ( |2.17| ) with Ew = W (for a Galilean boost). 

For X > 1 the two vectors pi and p2 always lie in the first or third quadrant, but for x < 1 

the vectors may lie in any quadrant. 

The restriction of both momenta pi and p2 to the first and third quadrant, which happens 

for X > 1, produces a curious singularity in the lab cross section. Under these conditions, 

the lab angle 61 will reach a maximum value less than 90° at a cm. angle 6* = 6*^.^^. At this 

point 



0. (2.20) 



The differential cross section in the cm. is always finite, but the lab cross section, defined 
by the transformation 



iMc.„..(^:nt)" 



da = de* |M,.^.(^;,j|^ = de, \M,,^{e:,,,)\' ^ de, ' ;^;^-^;gr' , (2.21) 



has a singularity at the critical point because of the vanishing of the Jacobian (|2.20|) . This 



singularity is analyzed in detail in Appendix |C|. In this paper we present cm. cross sections 
only, so we do not encounter this singularity. 
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FIG. 5: The symmetry momenta and the symmetry angles Og as functions of for five values of x: 
X = 2 (solid line), 1.8 (long dashes), 1.5 (medium dashes), 1.0 (dashes) and 0.5 (dots). The momenta rise 
with and the angles fall with Q^. For x ^ 2 the angle dg — and is not shown. 

For the special case when 6* = 7r/2 (where the relativistic elhpse touches the nonrela- 
tivistic circle) the magnitudes of pi and p2 are equal, and the RIA depends uniquely on the 
wave function at only one momentum point. We will refer to this as the symmetry point. 
Since the angle between pi and p2 is 90° in the nonrelativistic limit, this is referred to as 
perpendicular kinematics. Were there no final state interactions or interaction currents, the 
symmetry point would be an optimal place to measure the wave function. The symmetry 
point momenta ps = |pi| = |p2| and angles 9s = Oi = —62 are 



These are shown in Fig. ^ as a function of for several fixed values of x. The figure shows 
that if we wish to measure the wave function at large p^ (near one GeV) and at large x 
where pion production is not large, we must go to large (about 2 to 3 GeV^). 
We now turn to a discussion of the RIA. 

C. Matrix Element for the Relativistic Impulse Approximation 

The RIA approximation used in this paper is based on the simple pole diagrams shown in 
Fig. ^. We use a (standard but unfamiliar) notation in which matrix elements of an operator 
between two outgoing Dirac particles are written in the form 



where C = —i'-f^'j'^ is the Dirac charge conjugation matrix. This notation is very convenient 
because C transforms like an incoming v spinor (but is not to be interpreted as an antipar- 




(2.22) 




(2.23) 
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tide in this application), and therefore the most general operator O can be constructed from 
the standard 16 independent Dirac bilinear covariant operators. The matrix representation 
( p.23| ) is equivalent to a direct product representation 



«a(Pl, Ai) [OC]^. mJ(P2, A2) ^ Ma(pi, Ai)m^(p2, As) [OC]^ 



(2.24) 



but is more convenient for relativistic calculations. [Note that the RHS and LHS of this 
equation are identical as long as the Dirac indices are shown explicitly, but only the LHS can 
be turned into Eq. ( |^.23| ) by dropping explicit reference to the indices. Beware that the order 
of the momenta in ( |2.23| ) and ( |2.24| ) is opposite from that used in a previous reference 



where m(p2, A2) was multiplied from the left and m(pi, Ai) from the right; see Appendix ^] 
Including the isospin factor (H; | ~ | I 00^ = ^ the Feynman amplitudes for the RIA 
in any frame can be written 



(AiA2|J,(g)|A,) 



^i(Pi, Ai)jJ^^(pi,Pi 



V2Ar^ 

- M2(P2, A2)jf Hp2,P2 - g)rii^A,(Pl. ^) 



(2.25) 



where A, are the nucleon helicities and A^ the helicity of the deuteron, and is a normal- 
ization constant defined and discussed below. The subscript on the nucleon helicity spinor, 
Ui [suppressed in Eq. ( p.23| )], refers to whether it is particle 1 or particle 2, in the sense of 
Jacob and Wick [^ (see the discussion in Appendix 0). The nucleon current is 



m 



(2.26) 

and p — q off-shell. 



where p and p — q are nucleon four-momenta with p on-shell (p^ 
the superscript i = 1 (proton) or 2 (neutron), and the virtual photon has polarization vector 
Eg, where g = {0, s, a} with s and a the linear combination of photon helicities introduced 



in Eq. ( |2.10D [for more details, see Eq. ( [B14D ]. The relativistic deuteron wave function 



for a nucleon with momentum P — p off-shell and a nucleon with momentum p and helicity 
A on-shell (so that p^ = rn?) is defined to be ■ipx^x^{p, P), and is related to the normalized 
dnp vertex function F by 



[p-py 



N,T,{p,P)Cui{p,X)exAP). 



where the normalization constant 



2Md{2Txf 



-1/2 



(2.27) 



(2.28) 



is chosen to give the defined wave function a convenient normalization [see Eq. ( |A15| )]. 
The superscript (i) labels the choice of helicity convention (particle 1 or 2) for the on-shell 
particle, F^(p, P) is the normalized dnp deuteron vertex (with the off-shell particle on the 



left) first defined by Blankenbecler and Cook |3T|, and (P) is the deuteron polarization 



vector for a state with helicity and four-momentum P. We use the notation of Ref [|30| for 
F. Note that the normalization constant in ( |2.25| ) and ( |2.27| ) cancel; the Feynman amplitude 
depends only on the normalization of F and not on the convention used to normalize ip- Foi' 
further details, see Appendix |^. 
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D. The issue of gauge invariance 



The RIA is not gauge invariant by itself. This issue must be dealt with before we can 
proceed with the calculation. Here we discuss how this is done. 

Using the method of Ref. the RIA, together with final state interactions (FSI) 

and interaction currents (IntC), is part of a gauge invariant calculation. Once all of these 
pieces have been calculated and assembled, the result will be gauge invariant. Here we 
describe a convenient prescription that is (i) covariant, (ii) renders each of the individual 
contributions (RIA, FSI, and IntC) separately gauge invariant without altering their sum, 
and (iii) modifies each of the individual contributions as little as possible. The method was 
introduced in Ref. [Q, where it was also shown that the prescription guarantees that the 
RIA also gives the correct asymptotic result for deep inelastic scattering. 

If the individual contributions to the total current are denoted Jria, Jfsi, and JintC) then 
the prescription calls for each to be modified by the replacement 

Jx = Jx-^q-Jx (2.29) 

where X is any of the RIA, FSI, or IntC terms. Since the total current is gauge invariant, 
q ■ J tots! = 0, and 

Jtotal — "^total (2.30) 

SO the prescription does not modify the total current. Furthermore, since the photon helicity 
vectors are all orthogonal to g, ■ g = 0, 

e^- Jy, = e^- Jx (2.31) 

and the prescription has no effect on the contribution of each of the terms in the current. 
This prescription meets all three of the conditions listed above. 

Unfortunately, there is no uniquely correct way to modify the RIA so that it is gauge 
invariant. The choice proposed here is only one of many possibilities. 



E. Calculation of the structure functions 

The structure functions are obtained by squaring the matrix element ( |2.25| ) and summing 
over spins. There will be three terms: the two "diagonal" terms coming from the square 
of the proton term and the square of the neutron term, and the interference term. The 
diagonal terms can be calculated by expanding the density matrices 

^^Hp) = T.^'xliP, P) ® 41% P) (2.32) 

in terms of independent Dirac spin invariants, and then performing the sum over the off- 
shell particle degrees of freedom using Feynman trace techniques. The final result, given in 
Ref. [0, is a sum of squares of invariant functions and scalar products of four vectors, and 
is manifestly covariant. 
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In this paper we present an alternative method in which the structure functions are 
calculated by first expanding the off-shell nucleon in terms of on-shell nucleon degrees of 
freedom, and then computing the squares of the matrix elements. It is possible that this 
method will simplify the calculation of polarization observables planned for future work. 
Unfortunately, the results obtained using this method are not manifestly covariant (but 
they are, nevertheless, covariant), and, unless one is extremely careful, it is easy to make 
sign mistakes by dropping one of the many phases that arise when transforming helicity 
amplitudes. As a check of the results presented here, we have shown explicitly that our final 
analytical result for the diagonal term is identical to the result obtained in Ref. P^ . 



The discussion of this method begins by noting that the physical content of the matrix 
element (|2^ ) can be displayed by decomposing the off-shell nucleon into positive energy {u 



spinor) and negative energy (v spinor) states. For example, if we choose a four-momentum 
k = {El:, k}, then the states u{k, A) and 7^u(k, A) (which we use in place of the v spinors), 
with helicity A = ±|, are complete, and 

1 = J2 {^(k, A) u{k, A) + 7^ u{k, A) u{k, X)-f^} . (2.33) 

A 

(The definitions and normalization of the helicity states are discussed in Appendix 0.) 
It is important to realize that while this decomposition can be carried out in any frame 
using nucleon states with any on-shell four-momentum, the result may appear very different 
depending on the frame and the spinor states used to do the decomposition (even though 
the final numerical result will always be independent of these choices). 

In this subsection we record the results for the current Eq. ( ^.25] ) if the decomposition is 
made in terms of the states of the spectator nucleon (with four-momentum p2 for the proton 
term and four- momentum pi for the neutron term). The final result in the cm. frame [see 
Eq. Ol is 



(A,A,iJ,(,-)iA,)=^J-±i:i: 

V -LU/1 i\d xp A'jA^ 

|7^,(2A;) j1%^(p,^^*, go) 0fA|(P2) 41(^2 -^) d^^/lliu,) d^/l\u,) 

-^p(2A'2)jfA':,(P,^*,?o)0fM(pi)rfL^i,A,(^i) <a1(^i) 4^xM] , (2-34) 

where the d's are the rotation matrices, A = A'^ + Ag, Uj are the Wigner rotation angles 
resulting from the boost of the spectator nucleons with four-momentum pj from the lab to 
the cm. frame, and pj and 6j are the magnitudes of the on-shell spectator three-momenta 
and polar angles in the rest frame of the deuteron [i.e. the lab frame; recall Table 0]. The 
phase r]p{x) is 



and the matrix elements of the single nucleon current, in the cm. system, are 

Ja.a,,(P,^ ,go) - i ,..(p*,A,)jf(p:,p:-g*)7^«,(p-,A) if p = - ^''^^^ 
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TABLE III: Wave function combinations that enter the current. 





(n(p) + V2wip)) 


<^o (P) = 




_ m 




- \/3-i;s(p)| 




[V2u{p)-w{p)) 


4>I{v) = 


7s' 


, m 


(V2n(p)-u;(p)) 





with j = 1 or 2, but j ^ i- The deuteron matrix elements are defined in the deuteron rest 
frame using the expansion ( p.33| ), and are written 




1 
2 



(2.37) 



The 0's are combinations of the four scalar deuteron wave functions defined by Eq. ( ^.271 ). 
These are the S-state m, the D-state w, and the two P-state wave functions Vt and Vs, and 
expressions for the 0's are given in Table |ITT[ The final result ( ^.341 ) was obtained by boosting 
this result to the cm. frame, as shown in Appendix ^. 

This form of the Born term makes it easy to examine polarization observables, and gives 
a simple form for the unpolarized response tensors ( |2.11| ). Squaring the proton term [with 
(i) = (1)] and summing over spins (averaging over the initial deuteron polarization) gives 



R. 



'99 



^ /j Jt 

127rW \ ^ ^' 



m 



m^Md 
W 



-^99' 



{0O+(P2)}' + {<Pt{V2)Y\ + J~g, [{00 (P2)}' + {0r(P2)} 



+ + sin^2) [0^(P2)0O (P2) + 0^(P2)0r(P2)] [ (2-38) 



where the currents are 



7+ 

■^99' 


- Jig 


- „(!)+ 
Jig' 


+ J2g 


h ■(!)+ 
J2g' 












■^99' 


- Jig 


Jlg' 


+ J2g 


J2g' 












'^99' 


_ A^h 

- Jig 


Jlg' 


+ Jv 


Jlg 


_ „{1H 
J2g 


J2g' 


■ (l)^ 

J2g' 


- ,(1)- 

J2g 




•^99' 


_ A^h 

- Jig 


J2g' 


+ J2g' 


Jlg 


, „(iH 

+ J2g 


Jig' 


+ Jis' 


J2g ■ 


(2. 



Recall that g and (7' can be either 0, s, or a. The 12 individual current matrix elements are 



given in Table These exact expressions are easily evaluated, and the response functions 
determined from Table |T|. The square of the neutron term [with (i) = (2)] is obtained by 
replacing 1 <-> 2, and the result for the interference term is given in Appendix 0. 
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TABLE IV: Matrix elements of the current. All variables are in the cm. frame, 
and E — \J m? + p^, = psin^*, and p^ — pcos0*. For the proton current, 
substitute proton form factors for Fi and F2 and set the phase 5 = +. For the 
neutron current, j'^', substitute neutron form factors and set 5 = 

,r = . -L i., sin e- H ir = -J- F,p + 1F„„ cos 



3i 



0- s^-^{F,uo-2tF2E) 



Q 



_ sine* fj^ j^ yoE 



■ -a _ . 1 p gOP± 



32 



32 



mQ 

cose* 
V2 



m 



m 



FiE + \F2y0 



" = 5^- (FiE + \F2v^ 



F. The cross section in the quasielastic hmit 



We may use expression (|2.38| ) to look at the cross section at the quasielastic peak, where 
P2 = and a; ~ 1 (we assume here that Md = 2m). Near p2 = the minus components of 
the wave functions are both suppressed, and the leading contribution to the cross section 
comes only from the term proportional to 



+ 



IV? 



2 I 2 
u +w 



2 I 2 
U +W 



47rn(p2 



(2.40) 



where the momentum density is approximately normalized to 



n(p) d^p= — I (fp u'^[p) + ly^(p) 



(2.41) 



[the exact relativistic normalization is given in Eq. ( [A15| )j. At the quasielastic peak 6* = 
and the cm. momentum p that enters the current matrix elements given in Table |V|is fixed. 
From Eq. ( p.l7| ) and the condition x = 1 we obtain 



P 



W^qo 



2Ew 2E\yMd 



2 



(2.42) 



This gives Rtt = and Rlt = 0, and the following simple formula for the coincidence cross 
section 



m 



■GUQ') + tGI{Q' 



+ 2TGlj{Q^)tan^e/2 



(2.43) 
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This can be compared to the cross section for scattering from a free proton, which is 

In both of these formulae, r = (5^/(4m^), and G^; and Gm are the famihar electric and 
magnetic form factors, related to Fi and F2: 

Ge = F,- TF2 

Gm = F, + F2. (2.45) 

III. PREDICTIONS OF THE RELATIVISTIC IMPULSE APPROXIMATION 
A. How good is the RIA? 

The low energy Bernheim (1981) data 0] are shown in Fig. ^ This figure shows (i) that 
the relativistic effects are large and (ii) at these low missing momenta there is no evidence 
for relativistic effects of higher order in v/c, nor for model dependencies coming from the 
difference between the Argonne V18 and Model IIB wave functions. The importance of 
relativistic effects in the d{e, e'p)n reaction has already been emphasized in Ref. |]T^, and is 



confirmed in our calculations. For this reason we will dispense with further nonrelativistic 
calculations, and present only calculations with relativistic effects included. 

Figure and the accompanying Figs. ^ and ||, also show that the low energy cross section 
is reasonably well approximated by the (relativistic) RIA. The agreement is at the level of 
±50% for a drop in the cross section of four orders of magnitude. We take this as evidence 
that the RIA is sufficient for the kind of crude survey carried out in this paper. 



B. Covariant RIA predictions for high 

We surveyed the d{e, e'p)n reaction over the wide range of kinematical conditions sum- 
marized in Table [V|. We did calculations at 4 different four-momentum transfers = 0.5, 
1, 2, and 3 GeV^, and 6 different values of the Bjorken variable x = 0.5, 1.0, 1.25, 1.5, 1.8, 
and just below 2, near the highest value accessible in the d{e, e'p)n reaction at a given Q^. 
[The value of x=2 can be reached only in the elastic reaction, not in electrodisintegration.] 
In addition to and x, Table ^ gives the transferred energy u, the magnitude of the trans- 
ferred 3-momentum = |q^|, and the range in missing momentum, where p^*" corresponds 
to a value of 6* = 0°, and p™"^ corresponds to a value of 6* = 180°. In addition, we list the 
value of the final state np relative energy (in the cm. system), Enp, which is 



Enp = ^Jiv + MdY - ql- rup- mn, (3.1 
and the kinetic energy of the pn system in the Lab frame, 

.La, Mj , 

,mx 

The highest accessible value of x for a certain is given by 



T^:' = ^-2m + ^{^-l] . (3.2) 

F 2m 2m \i 
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FIG. 6: The Bernheim data at low missing momentum. Note that the relativistic effects (mostly from 
the current operator) are significant. [Three of the lowest p„i data points, possibly contaminated by brem- 
strahluung, have been omited from the figure.] 



which comes from the requirement that T^^'' > 0. In Table |V], we list the kinematic variables 
for X = Xmax the last line for each Q^. The closest x value to that is the highest one for 
which we present calculations later on; it is characterized by a kinetic energy in the Lab 
frame of roughly 10 MeV. 

Table [V| shows that the values of the transferred energy and transferred three-momentum 
are closest for low x. In nonrelativistic reduction schemes one often assumes u << g^, which 
is clearly not the case for low x. Note also that the np relative energy is highest for low x. 
These imply that relativistic effects should be very strong in this region. This is interesting 
as it allows for a description of final state interaction by Glauber theory . Although final 
state interactions (FSIs) are not considered in this present study, it is useful to know that 
they can be calculated reliably in this kinematic region. 

At high X, the np relative energy is very small, on the order of a few tens of MeV. 
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FIG. 7: The acceptance averaged relativistic RIA calculation compared to the Bernheim data. 

Under these kinematic conditions, the system is reminiscent of a bound system, and one 
might reahstically expect wave function physics to be important here, e.g. the presence of 
the relativistic P-waves. The calculation of final state interactions proceeds by including 
the lowest partial waves. So, we have reliable methods for the calculation of FSIs both at 
high X and at low x. In the region in between, the description of FSIs is more involved and 
accordingly more difficult. 

The case of x = 1 roughly corresponds to the quasi-free case. Strictly speaking, the 
quasi-free case corresponds to 



V 



2m 



-Eh 



(3.4) 



N 



where the binding energy E^, leads to a small deviation from x = 1. However, the binding 
energy for the deuteron is small and we will refer to x = 1 as the quasi-free case in the 
following discussion. 

C. Six approximations 

Six different theoretical approximations will be discussed in the following. The first three 
are based on the covariant spectator RIA presented in detail in this paper, and will be 
denoted "C-IIB", "C-IIB-noP", and "C-AV18". 



C-IIB is calculated using the covariant IIB deuteron wave function obtained from the 
successful IIB A^A^ interaction This wave function and the relativistic spectator 



model have been previously used to successfully explain the elastic deuteron form 
factors [^. The full wave function has 4 components: the familiar S- and D-states, 



and two small P-states of relativistic origin. 
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FIG. 8: The ratio of the Berhehii data to the acceptance averaged relativistic RIA calculation shown in 
Fig. [7]. The triangles are the low pm data set, and the circles the high pm set. 



• C-IIB-noP is calculated using covariant IIB S- and D-state wave functions, but set- 
ting the small relativistic P-state components to zero. 

• C-AV18 is calculated from the covariant spectator RIA formulae using the S- and D- 
state Argonne V18 deuteron wave functions (instead of the wave functions derived 
from the IIB one boson exchange model) and setting the P-state components to zero. 
This is included among the covariant models even though, strictly speaking, the wave 
function is not consistent with the covariant formalism. This model is very similar to 
the "covariant model" previously discussed by Arenhovel fl^ , but they used ordinary 
spin instead of helicity, and the Paris deuteron wave functions instead of the (very 
similar) AV18 wave functions. 

The next three calculations are not consistently covariant, but they do use relativistic current 
operators. They all use the nonrelativistic Argonne V18 wave function. The first of these is 



based on the work of Adam and Arenhovel |36 



AA-v/c uses a current operator that results from a f /c-expansion of the intrinsic 
current |3^, |3^. Matrix elements of this current are made frame independent by 
replacing the approximate noninvariant effective 3-momentum transfer derived in |3^, 
P7| by its invariant extension, defined to be: 



'ieff 



1 + 



+ (ed - Wf 



(3.5) 



16m^ 



where e^, W are energies of pn pair in its respective cm. frame in the initial and final 
states. This prescription is similar to one proposed long ago for the elastic scattering 
by J. Friar ||3^, but differs from calculations by Arenhovel et al ||17|. 
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FIG. 9: The differential cross section in the cm. system at = 0.5 GeV^ and for x — 0.5, 1, 1.25, 1.5, 1.8, 
and 1.89. Each panel shows the six calculations described in the text: C-IIB (long-dashed line), C-IIB- 
noP (solid hne), C-AV18 (dashed line), AA-v/c (widely dotted Hne), JD-fuU (dash-dotted hne), and JD-lst 
(closely dotted line). Note that the different approximations are hard to distinguish, as discussed in the 
text. 
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FIG. 12: The differential cross section at = 3 GeV^ and for various x. The meaning of the curves is the 
same as in Fig. ^. 
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TABLE V: Overview over the kinematics employed in the calculation of the differential cross section 
and asymmetry A^. 
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The next two versions are from Jeschonnek and Donnelly . 



JD-fulI uses a fully relativistic, positive energy current operator. This covariant 
current differs from the spectator one by certain off-mass-shell extensions studied in a 
recent paper by two of the authors P^|. 



JD-lst uses a current operator expanded to first order in the initial nucleon momen- 
tum, with all other terms retained fully. This approximate "first order" form should 
be closer to the covariant one than the traditional v/c current mentioned above, since 
an expansion is made only in terms of the moderate momenta of nucleons in the initial 
nucleus. 
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The relativistic one-nucleon current used here in the JD-full calculation has been recently 
employed by Donnelly in studies of (e, e'N) reactions. In these studies relativistic models 
appeared to be far more successful than nonrelativistic ones |]I7| , p9| . It is, however, a non- 
trivial task to extend them beyond RIA. 

Final state interaction and meson-exchange currents (MECs) have been so far included 
into realistic calculations mostly within approximate frameworks based on various expan- 
sions of the nuclear operators in terms of supposedly small momenta PUI EDI, ffll] . We do not 
intend to give an exhaustive survey of those techniques, neither do we dare to compete in 



completeness and consistency with recent elaborate calculations 



We only show in 



our figures the results obtained with the various one-nucleon currents introduced above. 

While a much more comprehensive study of relativistic effects, including relativistic ex- 
pansions of TT exchange currents and heavy meson exchange currents including boost terms, 
77rp and 7770; currents, and isobar contributions, was performed by Ritz et al [^] for lower 
energies, we focus on high energies. Here, high energies mean the GeV region, accessible by 
CEBAF at Jefferson Lab and even the new kinematic regime opening up with the planned 
12 GeV upgrade of CEBAF. Our C-IIB calculation given here is fully covariant, and part of 
a consistent treatment of the nuclear dynamics and the one-body current that does not rely 
on any kind of nonrelativistic expansion. 

None of our calculations is complete, the purpose is rather to explore various experimen- 
tally feasible kinematical regions to find those for which the complete microscopic calcula- 
tions and precise measurements would be worthwhile. Nevertheless, the variations between 
results obtained with versions of covariant currents and their approximations, as well as 
those between two covariant or two approximate formulations themselves, should provide 
some insight on the region of validity of the expansions and approximations used. 



D. Differential cross section 

The differential cross section ( |2.1| ) is given in Figs. |^12. The six panels in each figure 
all have the same scale, so the relative size of the cross section may be seen at a glance. 
Examination of these figures shows that the magnitude of the cross section depends strongly 
on Q^, x, and 9* . The bulk feature of the differential cross section consists of the two 
peaks at 9* = 0° and 9* = 180°. The first peak corresponds to the impulse approximation 
contribution, where the photon couples to the proton which is detected later on, and the 
second peak at 180° corresponds to the Born contribution, where the photon interacts with 
the neutron. For a; < 1.5 the two peaks are well separated because (cf. Fig. ^ the nucleons 
have very different momenta at ^* ~ and 180°. In this case one of the two RIA contributions 
[recall Eq. ( |2.1tj| ) and Fig. ^ is much larger than the other. 

However, if we wish to probe the deuteron wave function at high momentum, we will seek 
the region near 9* ~ 90°, where both nucleons have nearly the same momenta and only high 
momentum components of the wave function can contribute. In this region the cross section 
is very small (and FSI will be large), reflecting the small size of the deuteron wave function 
at large momenta. At large x, and in particular near x ~ 2, the two diagrams will always 
have large momenta (because the relative momentum of the final state is low), and the cross 
section shows no sharp forward or backward peak. Here high momentum components of the 
deuteron contribute over the whole angular range. 

Figures PHT^ seem to suggest that, for x > 0.5, the differential cross section is relatively 



insensitive to the model used. However, this is largely an artifact of the log scales used in the 



26 



1 .4 



1 .2 



ratio 



1.0 



0.8 



0.6 





x=0.5 










' ■ \ ^^^^^ 

■j, ^- 

\ \ / 
'^-^ 




•'/ 
■'/ 
•"/" 
■ ' / 
/ 

/ 


1 1 \ , .'' . . 





30 60 90 1 20 1 50 1 80 

e 



1 .4 



1 .2 



ratio 



1.0 



0.8 



0.6 







x=1.25 















30 60 90 1 20 1 50 1 80 



1.4 - 



1 .2 



1.0 



0.8 - 



0.6 - 



x=1.0 



7^ 



^ . 



7 



30 60 90 1 20 1 50 1 80 



1 .4 - 



1 .2 



x=i.5 : 



1.0 ^fLi - " — ^ 



0.8 - 



0.6 



30 60 90 1 20 1 50 1 80 



1 .4 



1.2 - 



x=1.8 



1.0 — ^^^^ 



0.8 - 

0.6 - 




30 60 90 1 20 1 50 1 80 

e 



1.4 - 



1 .2 



x=1.89 



1.0 



0.8 



0.6 



30 60 90 1 20 1 50 1 80 

e 



FIG. 13: Ratios of the differential cross section in the cm. system for the cases shown in Fig. ^ {Q^ =0.5 
GeV^). Here five of the calculations are divided by the C-IIB-noP calculation: C-IIB (solid line), C-AV18 
(long dashed hne), AA-w/c (short dashed line), JD-full (dotted hne), and JD-lst (dash-dotted line). 
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FIG. 14: Ratios of the differential cross section at = 1 GeV^ and for various x. The meaning of the 
curves is the same as in Fig. ^ The estimated errors for a JLab measurement using existing equipment, 
shown on panels for x = 0.5, 1.0, 1.5, and 1.95, are discussed in the text. 
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FIG. 15: Ratios of the differential cross section at = 2 GeV^ and for various x. The meaning of the 
curves is the same as in Fig. ^ The estimated errors for a JLab measurement using existing equipment, 
shown on panels for x = 0.5, 1.0, and 1.5 are discussed in the text. 
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FIG. 16: Ratios of the differential cross section at = 3 GeV^ and for various x. The meaning of the 
curves is the same as in Fig. ^ The estimated errors for a JLab measurement using existing equipment, 
shown on panels for x = 0.5, 1.0, 1.5, and 1.98, are discussed in the text. 
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FIG. 18: The at = 1 GeV^ and for various x. The meaning of the curves is the same as in Fig. 
The estimated errors for a JLab measurement using existing equipment, shown on panels for x — 0.5 
1.95, are discussed in the text. 
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FIG. 19: The at — 2 GeV^ and for various x. The meaning of the curves is the same as in Fig. |[ 
The estimated errors for a JLab measurement using existing equipment, shown on panels for x — 0.5 and 
1.97, are discussed in the text. 
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FIG. 20: The at = 3 GeV^ and for various x. The meaning of the curves is the same as in Fig. |[ 
The estimated errors for a JLab measurement using existing equipment, shown on panels for x — 0.5, 1, and 
1.98 are discussed in the text. 
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TABLE VI: Sensitivity to relativistic effects (R) compared to estimated measurement errors (M) 
for selected and x. 



(GeV)2 






1 








2 








3 




X 


0.5 


1.0 


1.5 


1.95 


0.5 


1.0 


1.5 


1.97 


0.5 


1.0 


1.5 


1.98 


da R (in %) 


50 


20 


10 


10 


100 


50 


20 


10 


im+ 


50 


50 


20 


da M (in %) 


1 


1 


1-3 


5 


5 


1-10 


1-10 




10+ 


1-20 


1-20 


50 


R (absolute) 


0.1 


<0.1 


<0.1 


<0.05 


0.2 


0.1 


0.1 


0.01 


0.5 


0.2 


0.1 


<0.05 


M (absolute) 


0.1-0.2 






0.1 


0.2-0.5 






0.1-0.2 


<0.1 






>0.5 



figures, and to show more clearly the relative size of the different calculations, Figs. |T3|-|1G 
show the ratio of each approximation to the C-IIB-noP calculation. The relative variation 
is largest for low x (a;=0.5), but is significant at all a;, varying from about ±10% to ±50%, 
particularly near 9* ^ 90°, depending on the values of and x. These variations are sum- 
marized in Table [V^ for the larger values of Q^. Since the cross sections vary by many orders 
of magnitude, this model dependence is not large enough to prevent these calculations from 
providing a useful estimate of the size of the cross section over a wide range of kinematics. 

For a few choices of kinematics, we have estimated the size of the experimental errors 
that can be expected from a measurement of this reaction at JLab using existing equip- 
ment. We find, for many kinematics, that the experimental errors would be small enough 
to distinguish between the different theoretical models shown in the figures. Our estimates 
of the experimental errors are shown in Figs. 1^16 and in Table 0. In all cases we have 
examined, except possibly at the largest values of x at the largest Q^, we could distinguish 
these models from one another. Of course, the final state interactions and exchange current 
contributions must be calculated before one has a complete picture of this process, but our 
results suggest that such a calculation is likely to be worthwhile. 

To estimate these errors we assumed the measurement would be carried out in Hall A 
with the hadron arm of the HRS spectrometer pair placed either in the direction of the 
momentum transfer vector (data points with solid circles) or to the left (triangles) or 
right (inverted triangles) of q^. Each spectrometer setting is able to measure a range of 
angles 9* , with the settings to the left and the right of q^ able to measure larger ^*'s than 
the setting along q^ (which samples angles near 6'* ~ 0). The errors grow as 9* gets close to 
the limit of the acceptance of the spectrometer, and this explains the large errors at certain 
angles shown in the figures (see for example, the case when x = 1.5 and = 3, where 
there are large errors at 9* 25°, 9* ~ 50°, and 9* ~ 80°). This is clearly an artifact of our 
crude estimates, and could be removed by repositioning the spectrometers. The statistics 
are based on running for one day at each setting under normal JLab operating conditions. 

These graphs show statistical errors only. The estimates of the statistical uncertainties 
were made by acceptance averaging and radiatively folding [Q] the PWBA model of 
Jeschonnek and Donnelly ||39|. An alternate 3-pole parameterization of the MMD nucleon 



form factors |^6| and the Argonne V18 NN interaction |^ were used. The simulations were 
done using a realistic acceptance model for the JLab-Hall A high resolution spectrometer 
pair. A maximum beam energy of 4 GeV (except 6 GeV was used for the = 3 GeV^, 
X = 0.5 case), beam current of 100 /xA on a 15 cm liquid deuterium target, and measurement 
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TABLE VII: Cuts used to arrive at the statistical 
uncertainties. 



(GeV)2 


X 




-cut 


x-cut 


1 


0.5 


0.8 


-1.2 


0.4-0.6 




1.0 


0.8 


-1.2 


0.8-1.2 




1.5 


0.8 


-1.2 


1.3-1.7 




1.95 


0.8 


-1.2 


1.9-2.0 


2 


0.5 


1.8 


-2.2 


0.4-0.6 




1.0 


1.8 


-2.2 


0.8-1.2 




1.5 


1.8 


-2.2 


1.3-1.7 




1.97 


1.8 


-2.2 


1.9-2.0 


3 


0.5 


2.6-3.4 


0.4-0.6 




1.0 


2.6-3.4 


0.8-1.2 




1.5 


2.6-3.4 


1.3-1.7 




1.98 


2.6-3.4 


1.95-2.0 



time per kinematic setting of 24 hours were assumed. The cuts shown in Table |V1J| were 
used to restrict the simulations to reasonable intervals around the desired kinematics. 



E. The Asymmetry 

Next, we present our results for the asymmetry A^j), which is closely related to the 
transverse-longitudinal response Rlt] see Eq. ( p.l4|) . The numerical results are shown 



in Figs. 1^-20. The asymmetry is zero for 6* = 0°, then becomes negative, with an ex- 
tremum around 6* = 60° for = 0.5 GeV^ and x = 0.5. The extremum shifts to smaller 
angles for increasing x and to larger angles for increasing Q^. Then, the asymmetry changes 
sign and exhibits another peak around 6* = 140° for = 0.5 GeV^ and x = 0.5. The 
positive peak shifts to lower angles both for increasing x and Q^. The appearance of the 
positive valued part of A^j) depends on the presence of the Born graph contribution, the 
impulse approximation alone would only lead to one negative peak. Accordingly, when both 
processes start to interfere, i.e. for the highest x values, the minimum tends to wash out, 
especially for the situations where the peak around 6* = 180° has vanished in the cross 
section. 

One can see at first glance that the asymmetry is less sensitive to the differences in the 
calculations than is the differential cross section, except near x = 0.5. Perhaps the most 
interesting feature of these calculations is the irregular shape of the asymmetry at a; = 0.5 
for = 2 and 3 GeV^. At = 2 GeV^ both versions of the JD and the C-AV18 calculation 
develop an extra dip. At = 3 GeV^, the results for x = 0.5 develop even more structure, 
with the AA-v/c calculation having the opposite sign near 6* near 0°, and the C-IIB-noP 
showing an extra peak around 9 = 90°. At = 3 and x = 0.5, A^ could give unique insight 
both into the effects of relativity and different wave functions. By contrast, at large x and 
large the asymmetry is very small and not measurable with sufficient accuracy. 

The uncertainties in A^ were generated by propagating the errors in the cross sections, 
where the latter included statistical errors folded in quadrature with an overall 5% systematic 



36 



uncertainty. Further, to simplify the procedure, the integrated yields for protons within 
the right and left hemispheres about the momentum transfer direction were assumed to 
correspond to = and = vr. Finally, the values of the cross sections, needed to 
propagate the errors for A^, were taken from the point (i.e. not acceptance averaged) 
values; the statistical uncertainties in the cross sections were, of course, determined using 
the full acceptance. 



IV. CONCLUSIONS 

In this paper we have estimated the d{e, e'p)n coincidence cross section using the rela- 
tivistic impulse approximation (RIA). Our calculations span the range 0.5 < < 3 GeV^ 
and X from 0.5 to just less than 2. In this kinematic region, we find that the results are 
sensitive to different approximate treatments of the single nucleon current, and conclude the 
following: 

• Using equipment already in existence at JLab, it is feasible to measure the unpolarized 
coincidence cross section over this entire kinematic range. The asymmetry can be 
measured at small x where it is large. 

• The coincidence cross section is sensitive to the theory over the entire kinematic range, 
and it appears that measurements can be done to an accuracy sufficient to distinguish 
a large variety of relativistic models from each other, except possibly when both x and 

are very large. 

• The asymmetry is less sensitive to the theory, except at the smallest value of x = 0.5 
where measurements can easily distinguish between different theoretical models. 

To complete this preliminary study, we must add relativistic final state interactions and 
interaction currents that are consistent with the RIA. This is certainly feasible at large 
values of x, where the low relative momentum in the final state makes it possible to use 
existing relativistic A^A^ interaction models. It may also be feasible at low x, where the 
large excitation energy deposited into the final state may justify the use of a relativistic 
generalization of the Glauber approximation. 
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APPENDIX A: THE DEUTERON WAVE FUNCTION 



1. Helicity spinors 



Following the conventions of Jacob and Wick the helicity spinors for a particle 
with four-momentum p = {E, p} are obtained from spinors with four-momentum {m, 0} 
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P=P1 



00 = 71+0, 



-P=P2 



FIG. 21: The definitions of tlie momenta and angles di and 02- 



and spins up or down in the z direction. The state is first boosted along the i-axis until 
its momentum is {E,0,0,p} and then is rotated it in the direction of p. The Lorentz 
transformation that does this is therefore 



Sip,Q=nip)BiQ = SiAip,Q) 
where B is the Dirac operator for a pure boost along the 2;-axis 



Cp/2 



with tanh^p = p/E, and TZ is the Dirac rotation operator 

n{p)=n{^,9,-<p: 



g- j Eg <l>/2^-i Sa e/2^i S3 4>/2 



(Al) 



(A2) 



(A3) 



that takes the momentum from along the i-axis into its final direction. Using this trans- 
formation in the xz plane (0 = 0), the helicity spinors for particle 1 are defined as in 



Refs. 24, 47 



u{p,\) = Mi(p, A) = S'(]3,Cp)m(0, A) = Mi(p, A,6'i) 

-^^^^^)-(p'^'°)-(2rsth^KJ^^^'^^ 

^;(-p,A) = t;i(p,A) = -(-l)5-^C^i*(p,-A)=t;i(p,A,^i) 



■Ryi^ei) t;i(p, A, 0) = 7^^(^l) 7Vwi(P, A, 0) 



cosh |Cp 



X.(^^i) (A4) 



where Ti-yifi) = Tl{0, 6, 0), 6i is the polar angle of the vector p = pi (see Fig. ^T]), and 

/ cos \ i—sm^6\ . , _s 



Note that the u{p, A) and v{—p, A) of Ref. |^ are identical to the ui(p, A) and f i(p, A) of 
Ref. [^. Following Jacob and Wick |2^, the helicity spinors for particle 2 are defined as 
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= TZyiO^) V2iv, A, 0) = nyi92) 7Vn2(p, A, 0) = " ci'hic!^' ) ^"'^^'^ ^^^^ 

where 6*2 = vr + 6^1 is the polar angle of the vector — p = p2 (see Fig. |2l|) and m(— p, A) 
and f (p, A) are identical to the M2(p, A) and f2(p, A) of Ref. |2^. Note that the angular 
conventions have < < vr and have n < 62 < 27i. Useful relations, valid on the two- 
component subspace, are 

ny{92)xM = 2Ax_.(^i) 

ny{n)xM = 2Ax_,W, (A7) 
and, on the full four-component space, Ui and U2 are related by 

M2(p, A, 62) = 2A TZyiTc) Mi(p, A, 61) . (A8) 
[Note that 2A = (— 1)^^'^.] This last formula is useful for the applications in this paper. 



2. Deuteron wave functions in the rest frame 

The deuteron wave function (|2.27| ) is manifestly covariant, and we use this feature to 
simplify the treatment. Applying Eqs. ( [A4|) and ( |A^) in the deuteron rest frame, the spinor 
for either particle 1 or 2 can be written in terms of the spinor with the momentum in the z 
direction 

Ui{p, A) = Ui{p, A, Oi) = TZyiOi) Ui{p, A, 0) . (A9) 

Note that this equation holds for particle 2 even though, in applications, we restrict 62 > vr. 
The active rotation of the deuteron helicity vector from an initial direction along the ;z-axis 
to an angle 6 with respect to the i-axis is given in terms of the spin 1 rotation matrices 

e.,W = E<\:,we„(o). (Aio) 

A' 

a 

Substituting (|A9| ) into ( p.27| ), working the operator TZ through the rest of the expression, 
and then using ( [AlOp to realign the deuteron helicity vector in the z direction by rotating 



it through angle —9, gives 

<i^(p,p) ^ ^«(p,^„) e (0) = TZyie,) EV^£(p>o) e (o)4?A,(^0 (ah) 



where we used d^xjx' ~ This argument works only in the deuteron rest frame 

where there is no total three- momentum to be rotated by TZy. 



In the original reference |30] the on-shell particle was taken to be particle 1 with four- 



momentum pi, and the wave function in the deuteron rest frame was expanded in terms 
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of (on-shell) particle 2 spinors p4|. [Also, be aware that the spinors used in Ref. pO| 



were 



quantized along the fixed z axis, in for v spinors the notation —s corresponded to spin 
projection +s in the z direction.] In the notation of Eq. ( |A6| ) this becomes 



A2 



(A12) 



Using the fact that ^2 = tt + ^i, so that 62 = n when = 0, and using Eq. ( All ) gives 



+t^2(pi,A2,^2)^AiA2,A'(Pl^' 



(A13) 



where the components i^'^^x2 x' (Pi^) have the relative momentum vector (same as the mo- 
mentum of particle 1 in the rest system) aligned along the +z direction. The ± components 
of the wave functions follow from the helicity spinors defined above and the expansions given 
in Ref. fS^: 



Kx^xSp) 
V'aiA^.a.Ip) 




^^(p) 

^t(p) 
' V2 



1(^2 

ia2 
72 



X., . (AM) 



Here u{p) and w{p) are the momentum space radial wave functions for the S and D states 
and Vt{p) and Vs{p) are triplet and singlet P state wave functions, which appear in a manner 
similar to the lower component wave functions for the Dirac equation. These wave functions 
are functions of the variable p = |p| and satisfy the normalization condition |^| 

p^dp {u\p) + w\p) + v^ip) + t;2(p)} = 1 . (A15) 

The wave functions ( |A14|) can be simplified by specifying the helicity states of the 
deuteron (in the rest frame). Since the deuteron is a particle 2 in the sense of Jacob and 
Wick, its helicity states are 



^^(±1,-2,0) ifAd = ± 



V2 



(0,0,1) 



ifArf = . 



Note, for future reference, that 



(A16) 



(A17) 



where 7^^ = lZy{Ti). If p is in the +z direction, it is not difficult to evaluate ( |A14| ), giving 

1 



'^AiA2,Ad(P^ 
'^AiA2,Ad(P^ 



»7r 

2Ai 



(^Ad.Aa-Ai /|Arf|(P) 



JAd, A2-A1 J|A^| 



i/|A.|(P) 



(A18) 
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where 



f^ip) = + ^/2w{p) /o-(p) = VSv^ip) 

/+(p) = V2u{p) - wip) /r(p) = VSvtip) . 
Combining the expressions (|A13|) and ( |A18|) gives 



(A19) 



A2 



-2Ait;2(Pl,A2,^2)/|A2-Ai|(Pl' 



(1) 

A2 — Ai, Arf 



(A20) 



One disadvantage of the expansion ( |A20| ) is that the four- momentum of particle 2, which 
is P2 = {Md — Ep-^,—pi} in the rest system, is not the same as the four- momentum of 
the spinor U2, which is p2 = {-Ep^,— pi}. This difference can lead to confusion, especially 
because the four-momentum p2 is not one of the four-momenta that naturally occurs in the 
problem. In this paper we avoid this confusion, exploit the freedom to expand the off-shell 
particle in terms spinors with any four-momentum, and choose the four-momentum of the 
on-shell particle (i.e. the spectator) for the expansion. The advantage of this choice is that 
four-momentum used to describe the off-shell particle is now one of the naturally occurring 
momenta in the problem. Also, we find that the formalism is simplified if we use 7^ u instead 
of V spinors to describe the negative energy states (as in Ref. |^). With this choice, we 
find the following expansion for the wave function 



/ 3 

V^Sa,(Pi,^) = \/^Ehl(Pl'^2,^l)0^,,+Ai|(Pl 



-2X2 7^ Mi(pi, A2, ^1) 0|A2+Ai| (Pl 



d 



(1) 



- Ao — Ai , Arf 



(A21) 



where and extra factor of -\/3 has been introduced for convenience. Projecting out the 
independent components, gives 



4(P) 
</'f(p) 



E. 



V3 



m 



En 



V3 



m 



00 (P) 

0r(p) 



4(-/onp)-/o-(p) 

^(-/i^(p) + /r(p) 



(A22) 



where the new wave functions were given in Table |T1[ 

Comparison of the expansions ( |A2CI| ) and ( [A21| ) underline the fact that the separation of 
the wave function into positive and negative energy parts is a matter of convention; only the 
total result is independent of this separation. 

By a similar argument, we expect the expansion for ip^'^^ to be 



^A2 A, 



:iP2,P) = \/^Eh(P2,Ai,02)0|A2+Ai|(P2 



-2Ai7^ M2(P2,Ai,^2)- 



'IA2+A 



i|(P2 



d 



S+A,A,(^2-vr) 



(A23) 
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One can prove that 



by finding the expansions for if)^'^'^ directly from those for 



Using (|A23|) and orthogonality relations for the spinors gives 

= M2(P2, Ai, 62) TxSP2, P) C M^(P2, A2, 62 



4A1A2 ni(p2, Ai, e^) rA,(p2, p) c n-' ui (p2, A2, e 

4AiA2( 



Ml(p2, Ai, 6*1) r_Ad(Pl, P) C Ml (P2, A2, 61) 



-1 



Mi(p2, Ai, e,) r„A,(pi, P)Cui (p2, A2, ^1) , (A24) 



where A = Ai + A2, Txj^{p,P) = r^(p, P) ^^^(P), Ry^n) p2 = pi, and we used the relations 
and (|A17| ). However, from (|A21D 



= {-l)^-'''u^{p„X2,0l)^^S,-x^PuP) 

= i-l f-^" Mi(pi, A2, ^^i) r_A,(pi, P) Cm[(pi, Ai, 0i) .(A25) 

Since this equation is symmetric under the interchange of Ai and A2, and pi = p2 in the cm. 
system, the two equations (|A24|) and ( |A25|) are equal, and 0+ = + . A similar argument 
holds of (f)~ . 

The wave functions used in this paper were obtained by solving the Spectator equation 
using a kernel adjusted to fit the A^A^ data below 350 MeV lab energy 



3. Boosting helicity spinors 

In order to boost the spectator equation wave functions it is necessary to have expressions 
for the pure boosts of the helicity spinors. Indeed, for the applications to elastic deuteron 
electromagnetic form factors and the response functions for deuteron electrodisintegration, 
it is only necessary to study the case where the boosts are made along the z-axis. Since the 
Dirac spinors for arbitrary momentum are defined in terms of a Lorentz transformation of 
the rest frame spinors, it is useful to define the four-momentum in the particle rest frame as 

p={m,0). (A26) 

Three basic Lorentz transformations are required in the following discussion. In the notation 
of Eq. (|A1|) , these are: 

P = S{pXp)p = 'n{p)}3{Cp)p 

p' = sip'Xp')p^np'mcp')p 

P = B{Qp (A27) 

The first Lorentz transformation connects the particle rest frame with the initial frame of 
the particle and is composed of a pure boost along the z-axis S(Cp) followed by a rotation 
into the direction of the initial particle direction TZ{p). The second Lorentz transformation 
connects the particle rest frame with the final frame of the particle and is also composed 
of a pure boost along the z-axis followed by a rotation. The third Lorentz transformation 
connects the initial and final particle frames with a pure boost along the ^-axis. 
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Now consider the boost of a hehcity spinor [either ui of Eq. ([A4|) or U2{p, A) of Eq. ( [A 61) ] 
along the 2;-axis 

= S{p', Cp') s-\p', B{Q Sip, Cp) «,(0, A) 

= Sip\ Cp') S {a-\p', Cp')BiCz)Aip, Cp)) «.(0, A) (A28) 

where S{p, (p) = S{A{p, (p)) is the representation of the Lorentz transformation A(p, (p) and 
the group composition property of the Lorentz transformations has been used in writing the 
final step. Since 

A-\p',Cp')B{QAip,Cp)P = P, (A29) 

this combination of Lorentz transformations must be equivalent to a rotation, and is referred 
to as the Wigner rotation. In this case, where the boost is along the 2;-axis and the momenta 
are in the xz plane, the Wigner rotation is a rotation about the y axis, denoted TZy{uJi). The 
boosted spinor can therefore be written 



B{Qu,{p,X) = ^(p',Cp')^s/K)wi(0,A) = ^(]5',Cp')E«i(0'^')4' 



(1/2), 



5:«,(p',A')4, 



(1/2), 



(A30) 



Since 7^ and C7° both commute with the boost, and since the boost is real, ( |A30|) also holds 
for 7^'Uj(p, A) and C -uf (p. A). 

In this paper B{C,z) will be chosen to be the boost from the cm. frame where P* = 
{Dq,0,0, —go} to the lab frame where P = {M^, 0,0,0}. This is accomplished by a boost 
in the positive z direction with tanh(^^ = qq/Dq. In the notation we have introduced, the 
basic equations are 



B{QCu^{v:A)=T.(^u^{v^A')d^^> 



(1/2), 



(A31) 



and the inverse relations 



B'\QCu"^{v^A)=T.Cu"^{v:A')d>l 

A' 



(1/2), 
A \ 



(A32) 



where the unstarred variables are in the deuteron rest frame (the lab frame) and the starred 
variables are in the cm. frame of the electrodisintegration process. 

The Wigner rotation angles can be computed from the standard relations 



B-\Q nyie,) BiCp.) = nyie*) s(c,)7^,(T^. 



(A33) 



where the upper sign holds for i = 1 and the lower for i = 2 [the change in sign is a 
consequence of the phase relation in Eq. (|A8|)], and 9^ = 9* and 6*2 = 6'* + vr [recall Table |; 
the form of ^2 is a consequence of the extra rotation by vr in Eq. (|A8|) ]. Writing the rotation 
and boost operators in closed form, Eq. (A33) can be written 

cosh - as sinh (cos ^9i - 2^2 sin i6'j) cosh ^Cp, - «3 sinh ^(p. 



(cos ^9* —iT^2 sin \9* 



cosh \C,p — 0:3 sinh \C,p {cos \uj + 1L2 sin \uj 



(A34) 
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where is a shorthand for either uji or —0^2, depending on the case under consideration. 
This operator relation can be separated into four independent equations relating u to the 
lab variables {pi, 9i] to the cm. variables {p,9*}. A convenient form of these equations is 



C-+ cos = ^2Md{E + m) cos ^{9* - u 



C__ cosi^i = J2Md{E-m) cosUe* + Q 



C++ sm\e, = ^2Md{E + m) sin 1(0* 
C+. sin {Oi = ^j2Md{E - m) sin ^{6* + u) (A35) 



where 



Cab = J{Do + Md){E, + am) + hJ {D^ - Md){E, - am) , (A36) 



with E = \Jm^ + jp- and Ei = ^m'^ + pf. The following identities, derived from these 
equations, are very useful: 

— sin 9* cos LJ — cos 9* sin uj 
m 

E ^ 

— sin 9* sin u + cos 9* cos uj 
m 

E ^ 

cos 9* sin UJ + sin 9* cos uj 

m 

E ^ 

— cos 9* cos UJ + sin 9* sin uj 
m 

—p cos uj 

psinoj = ——go sin (A37) 



E. . 
— sm 9i 

m 




cos 9i 




Do . . 
sm 9j 

Md 




DoEi 

cos 9i — 

Mdm 


qoPi 

Mdm 


Ei cos 9i — 

Md 


Do 

Md^' 


m . 

— — On sm 9i 





4. Boosting the deuteron wave function 

In order to calculate the response functions for deuteron electrodisintegration it is neces- 
sary boost the deuteron wave functions from the center of momentum frame of the final state 
proton-neutron pair to the rest (lab) frame of the deuteron (where the decomposition of the 
wave functions onto S, D, and P states has been defined). If the system is quantized such 
that the three-momentum transfer q lies along the 2;-axis, then the deuteron wave functions 
in the cm. must be boosted to the rest frame by a pure active boost B{(z) in the z direction 
with tanh^x = Qo/Dq, as defined in the previous section. 

The rest frame wave functions are obtained by applying the operator B = B{(z) to the 
wave functions ( p. 27] ), which gives 



= jlfpZ%2 IS-'} Nd T,M, P*) B-'} [BCuJipl X) 
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(1/2), 



(A38) 



where, in the next to last step, we used the boost properties of the hehcity spinors (|A31|) and 
fact that the propagator and Tx^{pi,P) are Lorentz scalars. Note that there is no Wigner 
rotation of the deuteron hehcity vector because the boost is in the same direction as its 
momentum (but the components of the vector do change). The wave function in the 
cm. frame can therefore be written in terms of the rest frame wave function 



(1/2), 



UJi 



(A39) 



Using the representations ( |A21| ) or ( |A23|) , and the boost formula (|A39|) and (|A30D , the wave 
function in the cm. frame becomes 




1A| 



B"' [n,(p„ a;.) </)+ |(p,) - 2a;. 7' a;.) 

A'A' 

J s 

[^i(P*.^)</'Jv|(Pi) -2Aj-7^Mi(p*,A)0^|(pi) 



A' A' A 



(A40) 



where A = A^ + A^-, 6'i = 9i and 62 = O2 — n, and the upper (lower) sign in the A index of 
is for i = 1(2). Note that the notation is mixed in the last equation : the momentum 
of the spinors is expressed in the cm. frame and the variables of the 0's and d^^^ are in the 
deuteron rest frame. 



APPENDIX B: THE HADRONIC MATRIX ELEMENT 



1. The plane wave matrix elements 



Using Eq. 



the current matrix element 



, in the cm. frame, becomes 



(A,A.|j.,(g)|A.) = ^ ' 



E 



Ml(Pi,Ai) jP(Pl,Pl -g*) M2(p;, A)0^|(p2) -2Ai7^M2(P2>A)0|^|(p2 



^2 



M2(P2, ^2)jf}{p*2,P*2 - q*) Ml(pt, A) 0^1 (Pl) - 2A2 7^ Mi(pt, A) 0|a|(Pi) 



iAr 



lA^ 



(Bl) 



where 



11(^2 -vr) d^y^i[{u:2) 41T(^2) 
.(^1) 4^a!K) <?(^i)- 



j(l/2) 



(1/2), 



(1) 

-A, A 



fB2) 
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and A and the Wigner rotation angles uji were defined above. 

The matrix elements of the single nucleon current operator, in the cm. system, are defined 
to be 



Ui 



.(P*, A.) - q*) M,(p*, A) if p = + 



Ui{v*Ad Jx! iP*i -(l*)l^Uj{p*,X) if p = - 



(B3) 



These are calculated in the next section. Using this notation, the current matrix elements 
( [Bl| ) can be written in the following compact notation 

AiA2|jA,(g)|A,) = (JAj^^,,(pt,p;,g*) 



3 1 



(1) 



7(1/2), 



;{l/2). 



(B4) 



where ?7p(x) is the phase defined in Eq. ( |2.35 ). The unpolarized cross section will be calcu- 
lated from this matrix element after the matrix elements of the nucleon current have been 
discussed. 



2. The single nucleon current 



In the spectator formalism used in this paper, the NN interaction kernel has a form factor 
h{p^) {h{m?) = 1) attached to each off-shell nucleon which enters of leaves the interaction. 
Alternatively, this form factor can be removed from the kernel and attached to the nucleon 
propagators, which then have the form 



m— 



(B5) 



Gauge invariance P2| will be insured if we introduce a reduced nucleon current jfi{p',p) 

rip',p)= h'Mp',p)h, (B6) 

where h = h{p'^) and h' = h{p''^), which satisfies the Ward-Takahashi identity using the 
dressed propagator 

qMp' ^P) = S-f\p) - ~Sf\p') ■ m 
A simple choice for the reduced current which satisfies this identity is PR R21 



3r{p\p) 



2m 

+Go(F3(g2) - 1) A4p')rA-{p) + GoA-{p')Y'A-{p) 



(B8) 



where Fi 2{Q'^) are the on-shell nucleon form factors, F^IQ"^) is a completely unknown form 
factor describing the off-shell structure of the nucleon (subject to the constraint that -^3(0) = 
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1), A_(j9) = (m— •f)j{2m)^ ;y/^ = _ q^^(^|q^ and Fq and Go are functions of and p'^ 
completely determined by the WT identity: 

1 — p'^ 1 — 



Since the final nucleon is always on-shell in the RIA approximation, p'^ = m? and the terms 
multiplied by A_(p') vanish, giving 

fnip'.p) = ^ [PriQ'W + F.iQ') ^ + ^} ■ (BIO) 

Furthermore, the terms proportional to g'^ vanish when the current is contracted with the 
photon helicity vectors. Hence, the current for use in the RIA reduces to the traditional 
current divided by h 

j''(p',p) = ^^{i^i(QV + ^2(Q')^^} {when p'' = m') . (Bll) 

Even when one of the particles is off-shell the only modification to the on-shell current which 
survives is the appearance of the factor of 1/h. 



The photon helicity states in the cm. frame, as defined in Ref. |21], are 



1 . 

^±1 = ^(0,Tl,-i,0) 
eo = ^(go,0,0,z/o) . (B12) 
Hence, the single nucleon current operator, defined in Eq. (|2.26|) and ( pil|) , is 

J^liPlM - Q*) = - , (B13) 

where the factor oi 1/h has been omitted (it is absorbed into the wave function). As 
discussed above and in Ref. the results are simplified if we consider the symmetric and 
anti-symmetric combinations of the transverse helicity amplitudes, which are found from 

Ss = ^(£i-£-i) = -^ (0,-1, 0,0) 

ea = \ (^1 + ^-i) = ^ (0, 0, -I, 0) . (B14) 

Using the explicit form of the spinors given in Eqs. (|A4D and (|A6|) we can show that the 
current matrix elements have the form 

(r. ft* n\ - \ ^A', -A Jlg^ + 2A' 5a',A jtg'' for {p, ^} = ( + , 0) , ( + , s) , ( - , tt) , . 

2A 5a', -A Jig + (^A'.A for {p,^} = (-,0),(-,s),(+,a) . 
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where g = {0, s, a} replaces the hehcity. If the proton (particle 1) matrix elements are 
calculated first, the neutron elements follow from 

Jyxt(P^O*,qo) = ^2(p,A',r + 7r)jf (p;,p;-g*)wi(p,A,r) 

= -4A'A Mi(p, A', r ) 7^42) (p;, p*^ _ q*) n-i ^^(p, a, r + vr) 

= -4A'A «i(p, A', e*)jf\pl,pl - q*) «2(P, A, e* + it) 

= -4A'AV(^7)jl}i+(p,r,-go), (B16) 

where the rotation by vr about the y axis has changed pi pi, q* ^ q* = (i/o,0,0, —go), 
£0 — — (— go, 0, 0, z/q) [so that the effect on the g = amplitude is to change the phase 
and to change go ~^ ~go], and therefore the phase rj'lg), arising from the rotation of the 
photon helicity vectors, is negative (for g = or s) or positive (for g = a). The change in 
the operator jjj^^ — > j^^^ corresponds to the replacement of the neutron form factors with 
proton form factors. The effect of the phase — 4A'A is an additional change in the sign of 
all _72 type elements. Combining all of these effects allows us to obtain jy x-^gip, go) from 

j?)^g{Pi^*,'lo) by changing proton to neutron form factors, go -go, ji(o,s) ^ -Ji(o,s), 
j2(o,s) j2(o,s), jia jia, and ^ -j2a- Thcsc phasc changes are recorded through the 
factor 6 shown in Table [TV] . 



The hadronic structure functions 



Substituting the form ( [B15|) for the current into the expressions (|B4| ) allows the sum over 
A to be carried out. If = or s the result is 

' ^ E E IvM^iK) [AT <I^U (P2) <-A. M + 2Ai AT (P2) d^lH (CO,)] y, 

^p(4A2AD [jiJ^^fAilPi) d^/Xi^r) + 2X,AT<f>\M(P^^ dx^xM)] ' (BIT) 



levr N, 



where 



3^2 

3^1 



4^AlK)rf-A,A,(^i)- 



(B18) 



lig 



a the phase 2Ai in the first square bracket multiplies j[T instead of j^T^ 



phase 2A2 in the second square bracket multiplies j[T instead of jig'*''- This different phase 
insures that there is no interference between and the other components of the current; 
Ja contributes only quadratically in the term J]. 

The unpolarized hadronic structure functions are now obtained by squaring the current 
( P17|) , summing over final hadron helicities and averaging over the initial deuteron helicity, 
and multiplying by the kinematic factors given in Eq. ( p.38|) . The structure functions i?*-^^^ 
are identically zero in the RIA. The others are proportional to 



T P 



A1A2 



(B19) 
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This generates three terms: the proton contribution (proportional to [j*-^-']^), the neutron 
contribution (proportional to [j*-^-*]^), and an interference term (proportional to 

j(l) X J (2)). 

The proton and neutron terms simplify easily. The sum over and A2 (for the proton 
term) or Ai (for the neutron term) collapses the sum over A'^ and A2 (associated with Jg and 
the sum over A" and A2 (associated with J^,) reducing the "diagonal" terms to 



X 



X'A' 



19 0^v|(Pi)+JS <^iAi(pj) 



+ 



(1/2), 



(B20) 



where i = 1 (for the proton) or 2 (for the neutron), j = 1 or 2 (but j ^ i), and we used 
: —AX[X d^^y\{uji) . Next, using the identities 



7(1/2) 



(1/2), 



7(1/2) , 



Y.<^i?i-^<^^xi-. 

X 



(1/2) 



rfL?,)(7r-2^,) 
0, 



26dlf)(2^,) 



(B21) 



and (2A)2 = 1, gives 



+ -^99' 



'P\A\{Vi) 



(B22) 



where the currents Jggi were given in Eq. (|2.39| ). Since these currents do not depend on the 
helicities, we may complete the sums using 



E Ci(p.) <i(p.) = 2 k(p,) 0^ (p,) + 0?(p,) 0f (p,) 



(B23) 



This gives the result reported in Eq. ( p.38|) . 

The interference term does not simplify as nicely. The sums over A^, Ai, and A2 can be 
carried out, but there are no delta functions to collapse the remaining four sums. The result 
is 



T 



12 



167r 



E '^A-A'(^2 



Ai At 



(B24) 



where A = A'^ + A2, A' = A" + A'2', and 



4" 



-2A 



d^'(^U^-) - ^KjS~ d^/^,{6.)] 0p^|(p2) 

^Kj[Td'xi!%is^)+j?J'^d'^yyU^4 *i(pi) 



(1)+ ^(1/2) 



+ 



+ 



",(2)" ^(1/2) 
Jig' « 



'2 

-5 ) - 2A" 7^'^" d^'/'^ 



'/'IA'iIpi 



(B25) 
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APPENDIX C: KINEMATIC SINGULARITY IN THE LAB CROSS SECTIONS 



In Sec. P^I B| it was shown that, if x > 1, the lab angle, 6i, reaches a maximum value 
in the first quadrant, leading to the condition ( |2.2CI|) . In this Appendix we show that this 
condition generates a true singularity in the lab cross section [suggested by Eq. (|2.21|) ], but 
that, because of the finite resolution of any detector, all observable cross sections remain 
finite. 

Mathematically, this singularity arises from a zero in the recoil factor TZ defined in 
Eq. ( p.8| ) . The denominator of Eq. ( |2.8| ) vanishes if 



Ew Pi = El cos 01 



(CI) 



where E\y = + v = \JW^ + is the energy of the final np pair in the lab frame. We 
will first show that the condition that this denominator vanish is identical to the condition 

To see this, it is convenient differentiate cos 6*1 with respect to 0* . Using Eq. (|2.17| ) gives 



d cos 9i d 

de* ~ de* 



Pi 

.Pi, 



pEw sin 9* pf 
2if 



p sin 9* 
p sin 9* 



E- 



w 



01 



q, p sin 9* + -^2p^ cos 9* sin 9* 



EwVi — Ql cos 9i 



Ew + ^pcosr 



{EwPi - q^ cos9iEi} 



where, in the last step we used 



Er = ]^Ew + ^Y>cos9* 



(C2) 



(C3) 



easily obtained form the same boost that gave Eq. ( p.l7|) . Equation shows that the 
two conditions ( p.20|) and (|C1|) are equivalent (except when sin^* = 0, when there is no 
singularity) . 

It is instructive to see how the cosine, 



COS^l 



p 



q^W + 2pEwz 



pf 2 + p2 sin^ 9* JqlW^ + M^pWEw Z + 4p2g2 ^2 ^_ 4p2^2 



(C4) 



of the proton lab angle, 9\^ varies with z = cos9*. For fixed x and Q^, depends only on z 
as given by Eq. (|C^). As an example. Fig. shows how varies with z for selected values 
of X when = 3. The value of z at which 2;^ is a minimum, denoted Zcrit can be computed 



from Eq. ( P4) . In the approximation that 



rrir 



m we obtain 



(iy2 + g2)(y^2_4^2) 



•^crit 



Wq, 



(C5) 



Moreover, solving Eq. ( PT| ) for z gives the same value, and hence we recover again the 
observation that the singularity occurs at the kinematic boundary of cos^i. 
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FIG. 22: The left panel shows — cos 6*1 as a function of 2 = cos 6** and the right panel shows \ — z^. In 
all cases — i GeV^, — 2 GeV, m — 1 GeV, and the lines show x = 0.5 (solid), x — 1.0 (very long 
dashes), x — 1.25 (long dashes), x — 1.5 (dashed), x = 1.8 (dotted), and x ~ 1.98 (dot-dashed). 



The values of Ei, pi and 2;^ at the critical point can be found using the additional 
constraint that follows from energy conservation. Energy conservation gives the following 
general formula for z^: 

_ ql+ml-ml + 2E^{Ma + u) - {Mg + vf _ 2EwE^ - 



where we have defined = W'^ — m\ + rn^. For fixed electron kinematics, Eqs. ( pi| ) and 
(PBD, taken together, give the proton energy Eq, proton momentum po and angle Zq at which 
the denominator (|Cl|) is singular: 



Of course, the value of Eq must be physical, i.e. the denominator has to be positive and 
Eq > rUp. The first condition leads to the constraint 

W'>ml-ml, (C8) 

which is always satisfied, while the second one requires 

g2 > 2z/(m„ - erf) - erf(2m„ - e^) ~ ^ , (C9) 

X 

where Md = nip + m„ — and in the last step we put m„ ~ nip ~ m, ~ 0, and 
u = Q'^/{2mx). This is the mathematical proof that the differential cross section is singular 
if and only if a; > 1. 

We conclude this discussion by showing that, even though the cross section is singular, 
the physical observables obtained from it are not. Because any measuring apparatus must 
necessarily have a finite resolution, all physical measurements must necessarily average the 
differential cross section ( |2.1|) over this finite resolution. At the kinematic boundary this 
average is 



dn'dE'dnj zoeJzo XdVl'dE'dVLi 
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This average will be finite only if the singularity is integrable, and this will now be shown 
explicitly. That it must be so follows from general physical considerations, and also from 
the behavior of the differential cross section in the cm. frame, where there is a smooth, 
non-singular behavior for all kinematical conditions. In the lab frame, for fixed electron 
kinematics, the dependent part of the integrand is contained in the kinematic factors 

p.'^{4n..] = ^-iS|_— ten...} (cii) 

At first glance, if we assume that the factor pi 7?. is an analytical function of z^, Eq. ( |(J1D 
suggests that the singularity will be a simple pole, which is not integrable. However, it is 
easy to show that the dependence of pi TZ on z^ is not analytic. Begin by using ( |C^ ) to 
rewrite Eq. ( P^ 

z^ = , a= , b = (C12) 

Pi Ql ^Ql 

and solve this equation for pi as a function of 2;^ (recalling that a and b depend only on the 
electron variables and are not functions of z^). The result is 



z^h + a^lP' + {zl - a'^)m'^ z^b + anip^'^ - Zq 
a? — 

where we used (|C7|) to express zq in terms of a and b 



Pi = ^ ^ = 5 — , (C13) 

— a-^ — zj^ 



la?m'i — 

zo = ^ . (C14) 



From ( |(J13| ) we obtain the following expression for the energy E 



f^pZr \ z'^-Zq + ba 

El = 2 • (C15) 

— z^ 



Substituting ( pi3| ) and ( pi5| ) into the recoil factor we obtain 



Pi ^ = F „ ^ = / , 2 • ^^^^^ 

J^w Pi - ^1 Ql TUp Jz\ - Zq 



From this it is clear that in the vicinity of the end point singularity at zq, the cross section 
Eq. ( |C10| ) behaves as Xj Jz\ — Zq and the singularity is integrable. 
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